Spectroscopic test of Bose-Einstein statistics for photons 
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Using Bose-Einstein-statistics-forbidden two-photon excitation in atomic barium, we have limited 
the rate ol statistics- violating transitions, as a fraction v of an equivalent statistics-allowed transition 
rate, to v < 4.0 x 10 -11 at the 90% confidence level. This is an improvement of more than three 
orders of magnitude over the best previous result. Additionally, hyperfine-interaction enabling of 
the forbidden transition has been observed, to our knowledge, for the first time. 



o 



43 
i 

a 

o 

-I— > 
O 
43 

(N 
> 



o 
o 

> 

X 



PACS numbers: 42.50.Xa, 82.50.Pt 

The Spin-Statistics Theorem (sst) is proved in the 
framework of rclativistic field theory using the assump- 
tions of causality and Lorentz invariance in 3+1 space- 
time dimensions. This letter describes an experimental 
test of Bose-Einstein (be) statistics and, consequently, 
the SST as it applies to photons interacting with atoms. 
As such, the experiment is a test of the assumptions, 
both explicit and subtle [Tj , of the SST and the assump- 
tions made in the quantization of the photon field. 

The experiment uses a selection rule [3] for atomic 
transitions that is closely related to the Landau- Yang 
theorem [H |5] in high-energy physics (hep). The selec- 
tion rule states that two collinear, equal-frequency, pho- 
tons cannot participate in any process that would re- 
quire them to be in a state of total angular momentum 
one. An example in hep is that the neutral spin-one Zq 
boson cannot decay to two photons, Zq -y* 77 (presently, 
the branching ratio for this process is limited to less than 
5.2 x 10 -5 [5]). For atoms, the selection rule means that 
two collinear equal-frequency photons cannot stimulate 
a transition between atomic states of total angular mo- 
mentum zero and one, J = -v> J' = 1 (Fig. [l]). The 
reason is the same in both cases: The photons would have 
to be in a state of total angular momentum J = 1, an 
exchange- anti symmetric state, in violation of be statis- 
tics and the SST (see Ref. [7] for a review of other ex- 
amples, other tests of particle statistics, and theoretical 
attempts to accommodate violations of the SST). We il- 
lustrate the selection rule, applied to this experiment, 
as follows: The rate W for E1E1 two-photon excitation 
from an atomic ground \g) to excited |e) state is [SI [H] 
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FIG. 1. (color online). Two amplitudes interfere in the exci- 
tation of the m'j — sublevel inaj = 0— >J' = 1 two-photon 
atomic transition. Similar diagrams can be drawn for excita- 



tion to the m'j — ±1 sublevels. 



is the amplitude for the ordered absorption of photons 
j and k; \g), |n), and |e) are the ground, intermediate, 
and excited states; r n and T are the intermediate- and 
excited-state natural widths; 61,2, ^1,2, and are the 
polarization, frequency, and field intensity of photon 1, 
2; LUki = (x)k—{x)i is the energy difference between states k 
and I expressed in frequency units; and £F is the electric- 
dipole operator. 

There are two paths to the final state, hence the two 

amplitudes an< ^ ^21 m ^Q- Bose statistics 

requires that the amplitudes add with a relative '+' (up- 
per) sign. We investigate the consequences of exchange 
antisymmetric photon states by adding the exchange am- 
plitudes with a relative ' — ' (lower) sign. Then, perform- 
ing the sum over magnetic sublevels of the largest con- 
tributing manifold, applying the Wigner-Eckart theorem, 
extracting the rank-1 irreducible component of the tran- 
sition operator, and assuming orthogonal light polariza- 
tions (which maximizes the rate |10j ) give 
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FIG. 2. (color online). Plot of /± 



on resonance (Sli + Q2 = 
The two 



u! eg ) as a function of one photon's frequency Oi. 
peaks correspond to the cases where the frequency of one of 
the photons matches the ground-to-intermediate transition. 
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FIG. 3. (color online). Partial energy-level diagram of bar- 
ium. Even (odd) parity states are in the left (right) column. 
Energies are in cm -1 . The close proximity of 6s6p 1 Pi to the 
energy midpoint (A = 18060 — 35934/2 = 93 cm - ) enhances 
the two-photon transition rate. 
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and where D nfl (D e n) is the reduced dipole matrix ele- 
ment of the ground- to intermediate- (intermediate- to 
final-) state transition. 

Equations Q and ^ show that when the two pho- 
tons have the same frequency (i.e., when fii = 572), 
then f + — 0, but /_ 7^ 0. Figure [2] shows this dis- 
tinction at f2 2 = fii = io eg /2: for exchange-symmetric 
photon states, the transition rate is identically zero; 
for -antisymmetric states, the rate is ocl/A 2 , where 
A = uj ng — uj eg /2 is the separation between the energy of 
the intermediate level oj ng , and the halfway point u> eg /2. 
This fact motivated the choice of atomic species in the 
experiment: in barium, A is relatively small (93 cm" 1 ), 
providing resonant enhancement to W- (FiG.[3]). 

Measurements of the transition rate made at laser fre- 
quencies near the halfway point probe the exchange sym- 
metry of photons. It is convenient to express the rates ^ 
in the case where the two laser frequencies are symmetri- 
cally detuned by a small amount 8 (where 8 <C A) about 
the energy midpoint, 8 = \uj eg /2 — = |w e9 /2 — Q 2 |. 
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expressed as functions of 8 are 



W- (8) 
W+(S) 



2-Dj n -Di g hI 2 

T/l 4 A 2 e2 c 2 



A 2 



A 2 



W- (6) . 



(6) 
(7) 



Because hamiltonians which treat multiple photons iden- 
tically cannot mix exchange-symmetric and exchange- 
antisymmetric subspaces of the photon state space [TT] , 
there is no interference between the two amplitudes in 
the total excitation rate. Thus the measured fluorescence 



signal S (5) on resonance is 

S (S) = 7 {W+ (5) + vW- (6)} 
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where 7 is the product of the number of atoms and the 
total detection efficiency, and v is a small dimcnsionless 
parameter that characterizes the degree of violation of 
be statistics. 

Measuring v is achieved by first calibrating the appa- 
ratus at (5 = 5 ca \. Expanding S (8) about 8 = 0, for 
v <C <5 2 al /A 2 -c 1, the calibration signal is 
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Measuring at the halfway point (5 
Sum = 5 (0) = ~fv W_ (0), so that 

= 5 Um 8^1 
5 cal A 2 ' 



0), the signal is 
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The previous version of the experiment [12] measured 
v < 10~ 7 at the 90% confidence level, using barium in a 
vapor cell and pulsed lasers to drive the 6s 2 -"So — > 5d6d 
3 S\ two-photon transition. As described below, the ex- 
periment's sensitivity was limited by the relatively large 
bandwidth of the lasers (« 3 GHz). This letter describes 
a new experiment that avoids these limitations by us- 
ing continuous-wave lasers, an in- vacuum power build-up 
cavity (pbc), and an atomic beam (Fig. [4]). 

A run begins with the lasers locked to pbc modes that 
are separated by 4 x c/2L = 2 GHz (Fig. [5) left). The 
cavity and lasers are then swept over LI, acquiring a cali- 
bration spectrum with <J ca i = 2ir x 1 GHz. The "Charyb- 
dis" laser is moved to the same longitudinal PBC mode 
as the "Scylla" laser. Then, both lasers are tuned back- 
and-forth over the limit point L2 many times, acquiring 
data with <5i; m < 27r x 2 MHz. During the «1 hour that a 
run can take, the lasers' intensities in the PBC may droop 
a few percent due to mechanical drift and dye consump- 
tion. The pbc transmission signals of both lasers are 
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FIG. 4. (color online). Experimental apparatus. An atomic 
beam of barium interacts with the light fields of two dye lasers 
(A = 557 nm) at the center of a tunable confocal optical 
power build-up cavity (pbc) of length L — 30 cm. The lasers 
are coupled into the cavity to excite only symmetric trans- 
verse modes, making the effective free spectral range (fsr) 
c/2L = 500 MHz. The cavity finesse is 340. The lasers drive 
the 6s 2 ^So — > 5d6d S D\ two-photon transition of barium, see 
Fig. [3] Fluorescence from 5d6d 3 Di to the 6s6p 3 P manifold 
(A = 420 nm) is collected through a light pipe and band- 
pass filters (not shown), and detected by a low-noise photo- 
multiplier tube in the photon counting regime. Each dye laser 
is locked to the PBC with an electro-optic modulator (eom) 
and rf-electronics, employing the Pound-Drever-Hall method 
|13| . Polarizing beam splitters and Faraday rotators isolate 
the lasers from each other and from their own cavity-rejected 
light. The lasers are also coupled into a single-mode fiber and 
interference is detected with a fast photodiode. A stabilized 
cavity, of FSR«50 MHz, records transmission peaks from the 
"Scylla" laser. 
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FIG. 5. (color online). Laser tuning path during run. The 
lasers are tuned in concert with the PBC, either separated by 
m—4: longitudinal cavity modes, or in the same mode. 



continuously recorded and used to correct for the droop. 
Scylla's transmission peaks through a stabilized reference 
cavity are recorded and used to calibrate the frequency 
axis and combine the data into a single spectrum, such 
as shown in Fig. [6j The stability of the atomic-beam 
density is monitored by measuring absorption on the 6s 2 
So — > 6s6p ^1 transition and found to vary negligibly. 

To extract an estimate of v from a run, a model line- 
shape Jzf(Oi) is constructed from the calibration spec- 
trum. The data in the limit region are then fit to the 
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FIG. 6. (color online). Processed data from one run. Fluo- 
rescence, in photons per second, from 5d6d Di to the 6s6p 



3 P manifold. {l[ = Hi 



,(138) 



/2 + <5 ca i, measured in ref- 



erence cavity fsr's (^50 MHz). In the region LI (Fig. IB), 
—4 < Cl[ < 11 , f2 2 — fii = 27T x 2 GHz, and a calibration 
spectrum is acquired containing both even and odd isotopes. 
Solid lines mark their approximate locations, known from a 
separate measurement to be presented elsewhere. Light-shifts 
are responsible for the broad line shapes. In the region L2, 
11 < il'i < 25, — f^i ~ 0, many traversals have been 
combined for improved signal-to-noise ratio. The spectrum of 
even isotopes (shown as dashed lines) is suppressed below de- 
tection by BE-statistics. In the spectrum of odd isotopes how- 
ever, the hyperfine interaction enables the transition. This is 
the first observation of this effect. 



model function aJzf '(f2i — O31/2) + b£l\ + c, where the 
parameters a, b, and c are determined by the fit, and 
fi3i/2 is the frequency difference between the calibration 
peak at LI and the limit point at L2. Without detailed 
knowledge of the light shifts, the exact frequency of the 
calibration peak and the limit point are uncertain. But 
all that is required is to know the frequency difference 
between them, and a simple measurement determines 
this: Tuning the lasers from LI to L3 (Fig. pi), and con- 
tinuously acquiring Scylla's reference-cavity peaks, two 
identical calibration spectra are acquired, one at LI and 
another at L3, separated by f^i = 41.5 reference-cavity 
peaks. Half this value is the frequency difference between 
the calibration peak and the limit point. 

The best-fit peak amplitude a is the ratio <5>ii m /<S C ai- 
Each run produces a best-fit estimate of v. The val- 
ues of v from twelve runs, the weighted average, and its 
standard error, are shown in Fig. [7j The best-fit value 
of v is (1.4 ±2.0) x 10~ n , limiting the relative rate to 
v < 4.0 x 10" 11 at the 90% confidence level. 

A number of sources of systematic error must be con- 
sidered. The Landau- Yang theorem is circumvented in 
the nonzero-nuclear-spin isotopes of barium, 137 Ba (11% 
natural abundance, I = 3/2) and 135 Ba, (6.6%, I = 3/2), 
by two manifestations of the hyperfine- interaction (hfi): 
HFl-induced splitting of the intermediate 6s6p : Pi state, 
and HFl-induced mixing of states with J' = 2 (primar- 
ily &sld 3 D 2 ) with the final 5d6d 3 D 1 state [HJ- The 
amplitudes for both processes are comparable, and the 
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FIG. 7. (color online), v from twelve runs. The error bars 
decrease over the lifetime of the experiment as improvements 
in excitation and detection efficiency were achieved. 



rate, relative to an allowed transition (i.e., expressed in 
the same way as the limit parameter u), is ~ 10 -8 . This 
rate is 10 3 times larger than the present level of sensitiv- 
ity, and the HFi-induced transitions are observed (e.g., in 
the right side of Fig. [6]). To our knowledge, this is the 
first observation of HFi-enabled two-photon transitions. 
Fortunately, the peaks of the HFi-induced transitions, al- 
though close to the limit region, do not completely ob- 
scure it. As is also discussed in Ref. [2], magnetic fields 
can create a false-positive signal in a way similar to the 
hfi. A relatively large field (~10 gauss) would be nec- 
essary to generate a false-positive signal at the present 
level of sensitivity. 

Another systematic effect, the one that limited the sen- 
sitivity in Ref. [T5] , is the non-zero spectral width of the 
lasers. Photons in opposite wings of the lasers' spectral 
profiles can together drive the nominally BE-statistics- 
forbidden transition at a relative rate ~ (r L /A) 2 ~ 
10~ 12 , where the laser line width is IY ~ 2tt x 3 MHz. 

When two counter-propagating light beams drive a 
two-photon transition, the Doppler effect due to the mo- 
tion of the atoms collinear with the light (i.e., the trans- 
verse motion in the atomic beam) shifts the frequencies, 
as seen by the atoms, of the two light fields in opposite 
directions. This leads to a nonzero allowed two-photon- 
transition amplitude. It is, however, largely suppressed 
in the case of a power-buildup cavity [10]. The residual 
effect has a relative rate T 2 Tp/A 2 ~ 2 x 10~ 15 , where 
T=10 -2 is the PBC mirror transmission, and rD=27rx 13 
MHz is the Doppler width of the two-photon transition. 

At the limit point, the amplitudes that destructively 
interfere (£^jk in Fig. [IJ proceed through magnetic sub- 
levels (or orthogonal superpositions thereof) of the in- 
termediate state. A disturbance of the degeneracy of 
these sublevels may alter the balance of the two am- 
plitudes, preventing perfect cancelation, and generate a 
false-positive signal in a way similar to the HFi-splitting 
mentioned above. Considering only the states in Fig. [T] 
light shifts produce a difference in the intermediate state 



Ei and Ei are the electric field intensities of the light. 
Notable is that both an imbalance in the light intensities 
and an imbalance in the dipole moments are necessary 
to make Si ^ 0. There are, besides, additional atomic 
states in the region around |e) that make smaller but 
significant contributions to Si. Roughly calculated, the 
relative-rate due to light shifts is (Si/A) 2 ~ 10~ 12 , an 
order of magnitude below the present level of sensitivity. 

The analysis of Eq. ([I]) and sequelae considered only 
E1E1 two-photon transitions. For the geometry of this 
experiment, where the light beams are collinear, the anal- 
ysis still holds true for all multipole combinations |10) . 
However, misalignment of the laser beams permits cer- 
tain higher-order multipole transitions [10]. The leading 
non-zero terms, E1M2 and E2M1 9 , roughly of the same 
magnitude, contribute a false-positive signal at a rate 
~ \1)M2/T)Ei\ 2 sm 4 9 ~ 10- 10 6» 4 ~ 10- 19 , where D M2 
and Dei are the magnetic-quadrupole and electric-dipolc 
reduced matrix elements of the transition, and 9 f=a 0.3° 
is a conservative upper bound on the misalignment angle. 

In conclusion, the reported experiment has improved 
the limit on possible Bose-Einstein statistics violation by 
photons by more than three orders of magnitude. Addi- 
tionally, we have observed hyperfine-interaction induced 
two-photon transitions. In principle, further improve- 
ment by an order of magnitude or more is possible with 
this technique by improving laser-frequency lock which 
will allow longer statistics accumulation. 
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